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Abstract. Proving the existence of an operator that connects non-perturbed states 
to perturbed states, an alternative derivation of the Dalgarno-Lewis method is given. 
To illustrate that the Dalgarno-Lewis method is an apt tool for algebraic Hamiltonians, 
the method is applied to one class of such systems, namely deep three-dimensional 
potentials with positive parity. 
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1. Introduction 

In perturbation theory the total Hamiltonian of the system under consideration is taken 
to be of the form 

H = H + Xh, (1) 

where it is assumed that the eigenstates and the eigenvalues of the Hamiltonian, Hq, 
are known: 

H \n) = e n \n) (2) 

and the parameter A is small. In this work we will assume that the eigenvalue spectrum 
of the operator H is non-degenerate. One is interested in finding the eigenvalue spectra 
of the operator H: 

#|* n (A)) = £„(A)|*„(A)). (3) 

The question we wish to address here is if there exists an operator S n (X) that 
connects each state in (T5]) and ([3]): 

\V n (\)) = S n (\)\n) (4) 

and describe the connection between such operators and the Dalgarno-Lewis method 
[TJ. Clearly one needs to impose 

Um5 n (A) = i. (5) 

We assume that the operator iS n (A) can be written in the form 

S n (X) = exp (AF n + X 2 G n + A 3 f n • • •) , (6) 
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where F n , G n and f n , • • • are operators to be determined. Inserting (j3J and (EJ) into (EJ), 
one finds that the condition 

([5„,(A)]- 1 ^(A)S n (A)) \n) = E n (X)\n) (7) 
needs to be satisfied. Expanding the eigenvalue, E n (\), in powers of A: 

E n (X) = e n + Ae« + \ 2 e^ + • ■ ■ , (8) 
one inserts into (JTj). After equating powers of A, one gets the following conditions: 

O(\):h-[F n ,H ] = tf (9) 

0{\ 2 ) : l -[F m [F n , H }} - [G n , H ] - [F n , h] = (10) 

0(X 3 ) : -[G n ,h] - [T n ,H ] + i[F ft , [G n ,H ]} + l -[G n: [F n ,H ]] 

+\[F n , [F n , h}} - \[F n , [F n , [F n , H }]} = (11) 
2 o 



and so on. The equalities in ([9]), (TlOl) . and (fTTj) are understood to be acting on the 
state |n), e.g. @ is a short-hand notation for (h — [F n , ^o])!^) — e n^\ n )- This equation 
suggests that we can take the operators F n to be the same for all n and the part of that 
operator which is diagonal in the basis of the eigenfunctions of H drops out from ([9]). 
Indeed multiplying (Q by the state (m| from left we find that 

(m\h\n) = (e n - e m ){m\F\n) + e$5 nm . (12) 

For n = m this equation yields the standard expression for the first order correction to 
the energy eigenvalues in non-degenerate perturbation theory. For n ^ m, f|T2|) gives 
non-diagonal elements of F. Without any loss of generality we will take the diagonal 
elements of this operator to be zero and write 

f=E(*)l* (13) 

We can now consider (Q an operator relationship provided that the right-hand side 
is replaced by a diagonal operator in basis of the eigenstates of H , the eigenvalues 
of which are e^. Following similar arguments we can also consider (TTU|) and fill I) as 
operator relationships provided that we replace the right-hand sides with appropriate 
diagonal operators. We will denote these operators e^\e^,e^ 3 \ • - •. We find it useful 
to spell out the diagonal and non-diagonal parts of the operator h: 

h = h D + A ND . (14) 
From OH]) we observe that hj) = and 

[F,# ] = W (15) 
Using f fT5|) . the operator equivalent of ffTUj) can be rewritten as 

- [G, Ho] - [F, h D ] - h ND ] = e^ 2 \ (16) 
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Inspection reveals that the first two commutators in (1161) yield non-diagonal operators, 
consequently one gets the standard second order result: 

~ 1 -(n\[F,h ND ]\n)=e^. (17) 

Dalgarno and Lewis developed a method for calculating second-order corrections to 
the perturbation expansion, which essentially uses the same commutator pQ. The non- 
diagonal elements of the operator G can be calculated by taking the expectation value 
of f|T6|) between states (m\ and |n) (m ^ n): 

( e (i) _ e (i)) A ii 
(m\G\n) = - _ m ' {m\F\n) - -—- (m\[F, h NB ]\n). (18) 

Again without any loss of generality we can take the diagonal part of the operator G to 
be zero. 

We now turn our attention to the third order correction. Using the Jacobi identity, 
after some algebra ( FlTl) can be rewritten as 

(-[G n , h D ] - [t n , H ] - ^[H , [F n , G n ]} + - [F%, H ] 

+^(F 2 h D + h D F 2 ) - PJiF^j \n) = e®\n). (19) 

Inspection reveals that the first four terms in ( Tl9l lead to only non-diagonal operators. 
Hence the third order correction to the energy eigenvalues is given by 

ei 3) = (n\l(F 2 h D + h D F 2 ) - F n hF n \n), (20) 

which is again the standard result. The non-diagonal part of ( TT9l can then be used 
to explicitly calculate the operator T. We omit that expression in this paper. We 
succeeded in showing that there exists an operator S that connects non-perturbed states 
to perturbed states in (TJJ and that this operator is the same for all states. Even though 
there is extensive literature on the Dalgarno-Lewis method [21 El IH [5], the existence 
of an operator that connects non-perturbed states to perturbed states, as well as the 
connection of this operator to the Dalgarno-Lewis method was not previously noted. 

Since the Dalgarno-Lewis method can be presented in terms of commutators, it is 
clearly very suitable for applications when Hamiltonians are formulated in an algebraic 
framework. We provide one such example in the next section. 



2. Rotational Bands in Deep Potentials 

If the Hamiltonians considered have an underlying algebraic structure, then calculation 
of the operators F, G, etc. are considerably simplified. As an application we consider 
a three dimensional Hamiltonian with a deep potential, which is a function of r 2 only. 
Such Hamiltonians arise in the study of alpha-clustering in nuclei. 
The low-lying states of the Hamiltonian 

H = ^ + V(r 2 ) (21) 
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can be evaluated by expanding the attractive potential in a power series in r 2 : 

V(v 2 ) = -V {1 - av 2 + (3r 4 + ■ • •) (22) 

Insight into the low-lying eigenstates of such a potential may be gained by taking the 
harmonic oscillator part 

H = + aV v 2 (23) 

as the unperturbed potential and 

h = -1/ /3r 4 (24) 

as the perturbation [6]. Perturbation theory can be formulated algebraically introducing 
the harmonic oscillator creation and annihilation operators 



1/2 / l \V2 

Xi -i — - p, (25) 



and 



\2hJ 

To utilize and algebraic formalism one introduces the operators K + , K_, K : 

^o = ^E( a la, + a4). (27) 

i 

K + = W*\*\> K- = {K+)\ (28) 
/ i 

which form the SU(1,1) Lie algebra with the commutation relations 

[K , K ± ] = ±K ± , [K+, K_] = -2K . (29) 
The Casimir operator of this algebra is 

C 2 = Kl- 1 - (K + K_ + K_K + ) . (30) 

For the realization given in (|27j) and f )28|) . this Casimir operator depends on the orbital 
angular momenta: 

C 2 = I L 2_A. ( 31) 
4 16 y ! 

The states are labeled by two numbers, k and m. The eigenvalues of the Casimir 

operator are 

C 2 \k,m) = k(k - l)\k,m), (32) 
which implies that 

k = l(e + l), (33) 



2 V 2 
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where I is the orbital angular momentum quantum number. The label m is the 
eigenvalue of the operator K : 

Ko\k, m) — m\k, m). (34) 

In the above equations m takes values of k, k + 1, k + 2, k + 3, • • •. One also has 

K±\k,m) = [(m± fc)(m=F fc± 1)] 1/2 |Jfe,m± 1). (35) 

It is easy to show that r 2 can be expressed as 

r 2 = —( K+ + K+2 K ). (36) 

fXU) 

Hence the perturbation Hamiltonian of (I24p is 

ft = -^-{K + + K_ + 2iT ) 2 = -^[6X 2 " 2C 2 + 2(K + if + K K + ) 

+2(K.K + K K^) + K\ + K 2 _\. (37) 

The first two terms in the above equation are diagonal in the unperturbed basis. It is 
then straightforward to verify that 



hp 



2(K_K + KqKJ) - 2(K + K + K K + ) + ^(K 2 _ - K 2 + ] 



2a^2afiV 

To calculate the second order correction to the energy eigenvalues, (|T7j) . we need to 
calculate the diagonal part of the commutator of F with Hnd- Defining the operators 

A+ = K + K + K K + , (39) 

and 

i_ = K^K + KqK- = (i + )t, (40) 

one sees, upon inspection, that there are two only diagonal contributions to [F,h ND ]. 
One comes from the commutator 

[Ah A.] = -16K 3 + 8C 2 K - 2K , (41) 
and the other one from the commutator 

\K\ , K 2 ] = -2(4K 3 - 4C 2 K + 2K ). (42) 
Putting it all together we can write down the energy eigenvalues of the Hamiltonian of 



2D as 



7mC 2 - 15m 3 - ^m) + ■■■ (43) 



2 J a 

where we defined the dimensionless expansion parameter 



s = n W» (44) 
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This Hamiltonian exhibits rotational bands (terms proportional to C<i ~ L 2 ) in its 
energy spectrum. The existence of such rotational bands for deep Gaussian potentials 
was proposed in [7] based on numerical calculations. We showed that application of 
the Dalgarno-Lewis method not only provides a convenient framework for calculating 
higher-order perturbative corrections, but also makes the existence of rotational bands 
manifest. 
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